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The subject of elastohydrodynamic lubrication is identified 
with situations in which elastic deformation plays a significant 
role in the hydrodynamic lubrication process. This chapter is 
concerned with the basic equations used in the analysis of hy- 
drodynamic lubrication. Basic elasticity theory is also consid- 
ered in this chapter, and elastohydrodynamic theory is developed 
in Chapter 7. 

Lubricants are usually Newtonian fluids, in which the rate 
of shear is linearly related to the shear stress. (A notable 
exception is grease, which behaves as a solid at shear stresses 
that are less than a threshold value and as a viscous fluid at 
higher stress levels.) The fluid is normally assumed to experi- 
ence laminar flow. Navier (1823) derived the equations of fluid 
motion for these conditions from molecular considerations and 
from the introduction of Newton's hypothesis for a viscous 
fluid. Stokes (1845) also derived the governing equations of 
motion for a viscous fluid in a slightly different form, and the 
basic equations are thus known as the Navier-Stokes equations of 


1 


motion. The Navier-Stokes equations will be derived as simply 
as possible. 

The study of hydrodynamic lubrication is, from a mathemati- 
cal standpoint, the application of a reduced form of these 
Navier-Stokes equations in association with the continuity equa- 
tion. The resulting differential equation was formulated by 
Reynolds (1886) in the wake of a classical experiment by Tower 
(1883) in which the existence of a thin fluid film was detected 
from measurements of pressures within the lubricant, Petrov 
(1883) had simultaneously recognized the existence of a coherent 
fluid film between the rotating shaft and stationary bearing in 
an extensive investigation of the friction of journal bearings. 
The Reynolds equation can be derived either from the Navier- 
Stokes and continuity equations or from first principles, provi- 
ded of course that the same basic assumptions are adopted in 
each case. Both methods will be used in deriving the Reynolds 
equation, and the assumptions inherent in reducing the Navier- 
Stokes equations will be specified. 

The Reynolds equation contains viscosity and density 
terms. These properties of the lubricant depend on temperature 
and pressure, and hence it is often necessary to couple the 
Reynolds equation with the energy equation. This chapter 
therefore deals with these lubricant properties and the energy 
equation. 
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The Reynolds equation also contains the film thickness as a 
parameter. The film thickness is a function of tne elastic be- 
havior of the bearing surface, and the governing elasticity 
equation is therefore presented in this chapter. The coupling 
of the Reynolds equation with the elasticity equation is consid- 
ered in Chapter 7, where the basic concepts of elastohydrody- 
namic lubrication theory are developed. 

5.1 Navier-Stokes Equations 

The Navier-Stokes equations can be derived from a consider- 
ation of the dynamic equilibrium of an element of fluid. It is 
necessary to consider 

(1) Surface forces 

(2) Body forces 

(3) Inertia 

5.1.1 Surface Forces 

Figure 5.1 shows the stresses on the surfaces of a fluid 
element in a viscous fluid. Across each of the three mutually 
perpendicular surfaces there are three stresses, yielding a 
total of nine stress components. Of the three stresses acting 
on a given surface the normal stress is denoted by o and the 
shear stress by t . To avoid overcrowding, the stresses on the 


3 


surface perpendicular to the 2 axis have been omitted. The 
first subscript on the shear stresses refers to the coordinate 
direction perpendicular to the plane in which the stress acts, 
and the second designates the coordinate direction in which the 
stress acts. The following five relationships should be noted 
in relation to surface stresses: 

(1) For equilibrium of the moments acting on the fluid ele- 
ment, the stresses must by syirwietric; that is, the subscripts on 
the shear stresses can be reversed in order. 

^xy ~ ^yx’ ^xz ~ ^zx’ ^yz ~ ^zy 

(2) The hydrostatic pressure p in the fluid is considered 

to be the average of the three normal stress components. 

°y + ^Z = -3p (5.2) 

The minus sign is used because hydrostatic pressures are com- 
pressive, whereas positive stresses are tensile. 

(3) The magnitude of the shear stresses depends on the rate 
at which the fluid is being distorted. For most fluids the de- 
pendence is of the form 



where 

n = a constant of proportionality known as the coefficient of 
absolute viscosity 
u = components of velocity vector 
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The terms in parentheses in equation (b.3) are a measure of the 
distortion of the fluid element. 

(4) The magnitude of the normal stresses can be written as 

Oi = - P + XaC + 2n (5.4) 

where 


_ _9u _9v 9w 

^ " 3x 9y 9z 


(5.5) 


and 

= a second coefficient of viscosity 

a 

The divergence of the velocity vector, the dilatation, 5 mea- 
sures the rate at which fluid is flowing out from each point; 
that is, it measures the expansion of the fluid. The factor 2 
is introduced in equation (5.4) as a consequence of the defini- 
tion of the coefficient of viscosity for a Newtonian fluid . A 
Newtonian fluid is one in which the shear stress is directly 
proportional to the rate of shear. This point is discussed in 
detail in texts on hydrodynamics where the relationship between 
shear and rate of strain components is considered in the 
direction of the principal axis (e.g.. Lamb, 1932). 

(5) By making use of equation (5.4) the following can be 
written 

Ojj + Oy + = - 3p + (3Xg + 2ri)C (5.6) 
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This equation shows that the average of the normal stresses 
differs from that previously describeo in equation (5.2). For 
these expressions to be compatible, the following must apply: 

= - I n (5.7) 

Tne forces due to the stress gradients must be added to the 
external body forces. Three stresses tend to move the element 
in the x direction. Thus the change of for example, 

across trie element and through a distance dy is ( 3 T^y 3 y)dy. 
This stress acts on the face of the fluid element with area 
dx dz and produces a force (3Tj^^/ay)dx dy dz. There are 
similar expressions for a and t . 

/\ X ii 

5.1.2 Body Forces 

The forces needed to accelerate an element of fluid may be 
supplied in part by an external force field, perhaps gravity, 
associated witli the whole body of tlie element. If the compo- 
nents of the external force field per unit mass are X , Y . 

3 3 

and Z^, these forces acting on an element are 

X^p dx dy dz, Y^p dx dy dz, Z^p dx dy dz (5.8) 
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5.1.3 Inertia 


The three components of acceleration of the fluid are the 
three total derivatives Du/Dt, Dv/Dt, and Dw/Dt. The signifi- 
cance of the total derivatives can be seen from the following. 
Consider the x component of velocity u. In general, 
u = f(x,y,z,t). The change in u that occurs in time dt is 
approximately 


In the limit as dt » 0, dx/dt = u, dy/dt = v, and dz/dt = w. 
Therefore, if equation (5.9) is divided throughout by dt, while 
making use of the above, we can write the total derivative for 
the u component as 


^ _ ju , ^ 

Dt 9t ^ 9x 


, 9u . 9u 
+ V -r — I- w — 
9y 9z 


(5.10) 


Similarly for the v and w components of velocity we can write 


Dv 

Dt 


9v 


9v 


9t 


9x 


... - . - . , 

— = — + u — + v-^+w — 


9v 


9z 


(5.11) 


Dw 

Dt 


(5.12) 


The total time derivative measures the change in velocity of one 
element of fluid as it moves about in space. The term a/3t is 
known as the local derivative since it gives the variation of 
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velocity with time at a fixed point. The last three terms are 
grouped together under the heading "convective differential." 

The mass of an element of fluid having dimensions dx, dy, 
and dz is p dx dy dz; thus the components of the resultant 
forces required to accelerate the element are 

P ^ 2 , p — dx dy dz, p dx dy dz (5,13) 


5.1.4 Equilibrium 


Having defined the body, surface, and inertia forces acting 
on a fluid element, we can now state the requirement for dynamic 
equilibrium mathematically. When the common factor dx dy dz 
is eliminated from each term and the resulting surface and body 
forces are related to the acceleration of an element of fluid, 
it is found that 


p 


Dt 


do 


9t 


px, + ^ + 

^ 9x dy 


9t_ 


xz 

dz 


(5.14) 


Also there are similar expressions for stress gradients and body 
forces that tend to move the element in the y and z direc- 
tions. 







+ 



(5.15) 


Dw 
^ Dt 


9t 


pZa + 


xz 


8t 


3x 


XZ 


3y 


3z 


(5.16) 
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The Navier-Stokes equations or the equations of motion can 
be obtained in terms of the velocity derivatives by substituting 
the results of equations (5.3) to (5.7) into equations (5.14) to 
(5.16). The resulting equations are 


p 


^ = 
Dt 


PXa 


_1R+ 2 A Jiu 
9x 3 9x I y9x 



,29 f du 
3 9x '^\9x 



+ 


9y 


3v , 3u 


3z 


3u , 3w 
9z 9x 


Dv „ i£ ^ 2 9 

p — = oY - -r^ 4* 

Dt ^ a By 


2 9 

~ flv _ iu\ 

+ 2 A 

3 9y 

V3y 

3 9y 


3v 

9y 


8w 

8z 


y (5.17) 


9 

f 9u 9v\ 

9 

r /9v 9w\ 

9x 

J^1^9y + 9x[ 

9z 

[_'’\9z + 9yj_ 


Dw _ 


+ 


2 _3_ 

/ 3w 9u\ 


/ 8w 8v\ 

3 3z 

^^3z 8xy 

3 3z 




The terms on the left side of these equations represent inertia 
effects, and those on the right side are the body forces, pres- 
sure gradient, and viscous terms in that order. Equations 
(5.17) are the most general form of the Navier-Stokes equations 
for a Newtonian fluid. 

For an incompressible fluid these equations reduce to 
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where 




D? ° "a - + 1 “ 


^ ^ ■ "’'a - I? + 


p = pz - -|^ + n w 
Dt a 9z 


V2 = ^ ^ ^ 

ax^ ay^ 3z^ 


In cylindrical polar coordinates with r ^ and 

c ’ c’ 

that X = r^ cos and y = sin 0^, equations (5.17) 
rewritten as 


/au 3u , V 3u , 3u v^\ 3n 


+ n y2 u 


r2 r2 

c c 


oi'57+ “^ + — ^ + »l^ + — I- PV -J--2E- 

^9t 3r^ Scp^. Sz r^J ” c 3t[i^ 


j. u2 V , 2 9u 

+ niv^ '' - — + :i — 

c ■'^c' 


r" 9cp 

C r T f 


Pl-^+ + 

■ c 


,9t ar„ ' acp^, ' " azy 


- If + n v2 „ 


(5.18) 


(5.19) 

such 
can be 


(5.20) 
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where 


v2=i + J:--i- + X^ + -iL 

9r2 Tj, 3r(. ^2 3^2 9^2 

c c 


(5.21) 


•"s* ®s» 

and 

<^s 

such 

sin i»s. 

and 

z = 

»'s 


cos Gg, the Navier-Stokes equations for an incompressible 
fluid are 


/ 9 u 5 u _v_ 3 u w 

Plat arg rg aOg + rg slxi 0g 3cpg 


8u 


= p^s + ^ - 


2u 2 3v 2v cot 6c 


3 w 


r2 r2 
s s 


-2 sin 0 3cp_ 

c* o o 1 a 


/ 8v 

pU,- + U 


9 v V 9 v 


w 3 v uv 

+ 


cot 0 


at atg tg 30g tg SIH 6 g Bcpg Tg Tg 


1 9 p o 2 8u 

S— -I. ^ nZ 1 — - 


2 cos 0 


s 8w 


" - rg aeg + n V V + ^2 80g ' ^2 gi^2g - ^2 gi„2e 3cPg 


/Sw , 9w , V 8w , 

) — + u H 1- ^ 

\9t 9r„ r 90« r sin 0 9cp r 


w ^ 4. ^ 


~ pZo 


1 / 2 


S rg sin 0g 3cpg 

2 cos 03 
r2 sin20 


n V w - 


w 


3 u 


r2 sin20 r2 sin 6 


(5.22) 
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where 


2 _ -L 


= 


ri 5=^3 


s 9r, 


r| sin 6s 


3 

86, 


sin 6 


s 


30 , 


+ 


sin^e acp2 
s s 


(5.23) 


5.2 Continuity Equation 

The Navier-Stokes equations contain four unknowns: u, v, 

w, and p. The viscosity and density can be written as func- 
tions of pressure and temperature. A fourth equation is sup- 
plied by the continuity equation. 

In continuous motion the increase in the mass of fluid 
within a fixed control volume in time dt must be equal to the 
excess of mass that flows in over the mass that flows out. Con- 
sider the mass flow across a fluid element such as that shown in 
Figure 5.1, The mass flowing into the element in the x direc- 
tion in time dt is (pu + [3(pu)/3xjdx}dy dz dt; thus the ex- 
cess of mass flow into the element in the x direction is 
-[ 3{ pu) / 3x jdx dy dz dt. Similarly, when the directions y and 
z are considered, the total excess of mass flow into the ele- 
ment in time dt is 

9 8 

fx Ty Tz 


12 


Now the original mass inside the element is p dx dy dz 
and the increase in mass inside the element in time dt is 
(3p/3t)dx dy dz dt. Thus 

(pu) (pv) + -^ (pw) = 0 (5.24) 

This then is the continuity equation in Cartesian coordinates. 

If the flow is steady, 3p/3t = 0. Also if the fluid is incom- 
pressible, p is constant and the continuity equation reduces to 


9u 8v 3w 
3x ^ 3y ^ 3z "" ® 


(b.25) 


5.3 Reynolds Equation 

The differential equation governing the pressure distribu- 
tion in a fluid-film bearing is known as the Reynolds equation. 
As pointed out earlier, this equation was first derived in a 
remarkable paper by Reynolds (1886). Reynolds' classical paper 
contained not only the basic differential equation of fluid-film 
lubrication, but also a direct comparison between his theoreti- 
cal predictions and the experimental results obtained by Tower 
(1883). Reynolds, however, restricted his analysis to an incom- 
pressible fluid. This is an unnecessary restriction, and 
Harrison (1913) included effects of compressibility and dynamic 
loading in a later analysis. In this section a generalized 
Reynolds equation is derived from the Navi er-Stokes and contin- 
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uity equations after assumptions that are normally valid for 
fluid-film luorication have been introduced. The Reynolds equa- 
tion is also derived directly from the principle of mass conser- 
vation and simple expressions for Couette and Poiseuille flow. 
The various terms in the Reynolds equation are discussed, and 
various forms of the Reynolds equation as they apply to particu- 
lar cases of bearing operation are presented at the eno of this 
chapter. 

5.3.1 Derivation of Reynolds Equation From the Navi er-Stokes 
and Continuity Equations 

The basic equation of fluid-film lubrication, the Reynolds 
equation, can be derived from the reduced form of the Navier- 
Stokes and continuity equations. Recall that a Newtonian fluid 
was assumed in the derivation of the Navi er-Stokes equations. A 
minimum number of restrictive assumptions are introduced that 
enable a general form of the Reynolds equation to be estab- 
lished. The assumptions are 

(1) Inertia and body force terms are negl igible compared 
with the pressure and viscous terms. 

(2) There is a negligible variation of pressure across the 
fluid film ( 3p/8z = 0) . 
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(3) Owing to trie geometry of ttie fluid film the derivatives 
of u and v with respect to z are much larger than the 
other derivatives of velocity components. 

The last two assumptions arise from the geometrical form of 
the lubricating film. If h and £ represent typical dimen- 
sions of film thickness and bearing length, respectively, it is 
found that the ratio h/£ rarely exceeds i0“^, is often as 
small as 10"^, and is most coir^rionly about 10“^. The lubri- 
cating film can thus be seen as a long, thin channel of small 
taper in all realistic bearing configurations. 

The application of these assumptions allows the havier- 
Stokes equations expressed in equations (5.17) to be reduced to 



(5.2b) 


Since the pressure has been assumed to be a function of x 
and y only, these equations can be integrated directly to give 
general expressions for the velocity gradients: 


3u _ z 
8z n 9x 



9v z 9p ^ 

8z ~ n ■ 3y n 


(5.27) 


where A and C are constants. 
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Now the viscosity of the lubricant may change considerably 
across the thin film (z direction) as a result of temperature 
variations that arise in some bearing problems. In this case, 
progress toward a satisfactory Reynolds equation is considerably 
complicated. 

An approach that is satisfactory in the majority of fluid- 
film applications is to treat n as the average value of the 
viscosity across the film. Note that this does not restrict the 
variation of viscosity in the x and y directions. This ap- 
proach is pursued here, but a more general derivation of the 
Reynolds equation that takes account of fluid property varia- 
tions along and across the film and that discusses in detail the 
assumptions made earlier has been presented by Oowson (1962). 

With n representing an average value of viscosity across 
the film, the velocity components become 


U 


= ^ 4. 

2n 9x 




A^ + 
n 


2^ 9p 

~ 2n 9y 


Z 


J 


(5.20) 


If we assume zero slip at the fluid-solid interface, the 
boundary values for velocity are 


z 


z 


- 0, u = u^, V = vjj 
= h, u = U3, V = Vg 






(5.29) 
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The subscripts a and b refer to conditions'on the upper 
(curved) and lower (plane) surfaces, respectively. Therefore 
u,, V,, and w, refer to the velocity components in the 

da d 

X, y, and z directions of the upper surface and Uj^, Vj^, 
and refer to the velocity components of the lower surface 
in the same directions. 

With the boundary values given in equation (5.29), the 
velocity gradients and components can be written as 


U 


^/ 2z - h\ 3p _ ub - 

3 z 2n / 3x h 

iX ^ /2z ~ h\ 3P ^ ~ 

3z y 2n / 3y h 



ua 











(5.30) 


(5.31) 


With these expressions for the velocity gradients and tne veloc- 
ity components we can now derive expressions for the surface 
stresses and the volume flow rate. 


Surface Stresses 

The viscous shear stresses acting on the solids as defined 
in equation (5.3) are (n 8u/32)^_g and -(n 3'^/3^)2=h' 

These expressions can be evaluated from equation (5.30) as 
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n iu\ _ h lE _ ~ ^ 




- 4 _ h 9p . - ^a) 

' ' 8z ^ '' ^ ^ 


z=h 


2 9x 


(5.32) 




Volume Flow Rate 


The volume rates of flow per unit width in the x and y 
directions are defined as 


q 


X 



and q 

y 



(5.33) 


Substituting equation (5.31) into equation (5.33) while Integra- 
ting gives 


q 


X 


qy 


JiL 

12n 9x 





9p 

12n 9y 



(5.34) 


The first term on the right side of these equations represents 
the well-known Poiseuille (or pressure) flow, and the second 
term represents the Couette (or velocity) flow. 

Returning to equation (5.31), the Reynolds equation is 
formed by introducing these expressions into the continuity 
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equation" (5.24) . Before doing so, however, it is convenient to 
express the continuity equation in integral form. 


f [ 1 ^ 


dz = 0 


Now a general rule of integration is that 


0 


3x 



f (x,y ,z)dz 


f (x,y,h) 


Ih 

9x 


Hence, if the density p is assumed to be the mean density 
across the film, the integrated continuity equation becomes 



■ »(“a - V ■ 0 (5-35) 

The integrals in this expression represent the volume rate 
of flow per unit width (q^ and q^) described in equation 
(5.34). When these flow-rate expressions (which are derived 
from the Navier-Stokes equations) are introduced into the inte- 
grated continuity equation, we obtain 


i ph^ 3p\ , 3 

( fihilpi, _L 

P(u^ + Ug)h“' 

. J_ 

"p(Vb + v^)h 

_ IL + 

\ 12n 3x/ 3y 


_ 2 _ 

3y 

__ 2 J 


+ P(w^ 




)-PUa 


ih _ pv Ih + h 
3x a 3y 



(5.36) 
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Tne first^wo terms are Poiseuille terms and describe the net 
flow rates due to pressure gradients within the lubricated area; 
the third and fourth terms are the Couette terms and describe 
the net entraining flow rates due to surface velocities. Tne 
fifth, sixth, and seventh terms describe the net flow rates due 
to a squeezing motion, and the last term describes the net flow 
rates due to local compression. The last four terms can be com- 
bined and written as 3(ph)/8t. The generalized Reynolds equa- 
tion then becomes 

J_ fph! ilV J_ fphl 9 , 8 9(ph 

Sx \12n 3x/ 3y Vl2n 8yy 9x |_ 2 J 3y [_ 2 J'*' 3t 

(5.37) 

Equations (5.36) and (5.37) represent the general form of the 
Reynolds equation within the range of assumptions discussed. 

5.3.2 Direct Derivation of the Reynolds Equation From the Laws 
of Viscous Flow and the Principle of Mass Conservation 

The Reynolds equation can be derived directly by consider- 
ing a control column fixed in space and extending across the 
film. Before discussing tne control volume, however, we note 
that three additional assumptions were introduced into the deri- 
vation outlined in the previous section. These additional as- 
sumptions were 
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(1) The flow is laminar. Turbulent flow is sometimes en- 
countered in high-speed, fluid-film bearings, and a modified 
Reynolds equation is required for this condition. 

(2) There is no slip between the fluid and bearing solids 
at common ooundaries. 

(3) The lubricant properties (viscosity and density) are 
constant across the thickness of the film. 

If these three assumptions, along with those given at the begin- 
ning of the previous section, are used, the two approaches to 
the derivation of the Reynolds equation will be compatible. 

Consider the rate of mass flow through a rectangular sec- 
tion control volume of sides ax, ay fixed in the coordinate 
system and extending across the lubricant film between the two 
bearing surfaces as shown in Figure 5.2. Note that one bearing 
surface is represented by the plane z = 0 ana the other by a 
curved surface such that the film thickness at any instant is a 

function of x and y only. This is exactly the coordinate 

system used in the previous derivation of the Reynolds equation. 

The mass of lubricant in the control volume at any instant 
is ph AX Ay. The rate of change within the control space 
arises from the change in both the density and the height of the 

column and is given .by [( 3( ph)/9t) ]ax Ay. The difference be- 

tween the rate of mass flowing into the control space and the 
rate of mass leaving the control volume is -(aq^/8x)Ax Ay in 
the X direction and -(9qy/3y)Ax Ay in the y direction. 
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The principle of mass conservation demands that the rate at 
which mass is accumulating in the control space (3(ph)/3t) must 
be equal to the difference between the rates at which mass en- 
ters and the rate at which mass leaves. Therefore 

- IT - Tf ■ "»') 

But 


and 



9h ^ ^ 9p 
P ^ 9^ 


_ 9 _ 

9t 


(ph) = p 



9h 

9x 



+ h 


9t 


(5.39) 


By making use of equations (5.39) and (5.34), equation 
(5.38) can be rewritten as 


_9_ /phi i£\ 9 /phi 92.\ 

9x yl2n 9x/ 9y yl2ri 9yy 


9 

Ph(Uj, + u^) 

9 

Ph(vj, + v^) 

3x 

L . 2 J 

9y 

2 _ 


+ p("a-»b--a|S-v,f) + h|£ (5.40) 


This is a generalized Reynolds equation for the assumptions dis- 
cussed and is exactly the same as that derived in the previous 
section (see equation (5.36)). The equation can of course be 
reduced to 

9 / ph^ 9p \ , _9_/ Ph^ _9£\ _ 

9x Vl2n 9xy 9yVl2n 9y/ 9x 



1 3 


2 

9y 

2 




to be consistent with equation (5.37). 


(5.41) 
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5.3.3 Standard Reduced Forms 


The specification of boundary velocities u^^, Vj^, and 
for the plane surface does not present any difficulty. In gen- 
eral the surface velocities of the upper bearing components need 
to be related to the translation of the component parallel to 
the coordinate axes and the rotation of the component about its 
own center. In subsequent chapters we shall be concerned with 
boundary velocities confined to the following values: 




U 

a 





= 0 


w 

a 


u 

a 


111 

8x 


(5.42) 


where u^ and u are the surface velocities of the two 
b a 

bearing components in tne x direction. Note that the coordi- 
nate X was associated with the semiminor axis of elliptical 
contacts in Section 2.2.4. If the surface velocities do not lie 
in the direction of the semiminor axis, the coordinate x 
should always be taken as the direction of surface motion, even 
if k is then less than unity. 

For steady-state conditions the Reynolds equation (as 

stated in equation (5.40)) becomes 

J_ ('£hl l£V+ _L f£hl = 12u (5.43) 
9x V n 8x/ 9y V h 3y/ 9x 


Where 


^a ^ 


u 


2 


(5.44) 



Again it shoulo be emphasized that equation (b.43) allows 


for variation of the viscosity and density in the x and y 
directions. If, however, the fluid properties do not vary sig- 
nificantly, the corresponding Reynolds equation is 

I?) + I?) = 12un (5.45) 


Equation (5.43) not only allows the fluid properties to 
vary in the x and y directions, but also permits the bearing 
surfaces to be of finite length in the y direction. Side 
leakage, or flow in the y direction, is associated with the 
second term in equation (5.43). If the pressure in the lubri- 
cant has to be considered as a function of x and y, the solu- 
tion of equation (5.43) can rarely be achieved analytically. 
Approximate numerical solutions are therefore sought as outlined 
in Chapter 7, where the elastohydrodynamic theory for elliptical 
contacts is developed. 

In many conventional lubrication problems side leakage can 
be neglected, and this often leads to analytical solutions. If 
side leakage is neglected, equation (5.43) reduces to 

(5.46) 


This equation can be integrated with respect to x to yield the 
familiar integrated form of the Reynolds equation 

ph - (ph). 


^ = I2nu 
dx 


'm 


Ph- 


is. 47) 
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where the subscript m refers to the conditions at points 
where dp/dx = 0, such as the point of maximum pressure. If the 
lubricant is incompressible, equation (5.47) reduces further to 


^ = 
dx 



(5.48) 


5.4 Viscosity 

The viscosity of a fluid may be associated with its resis- 
tance to flow, with the resistance arising from i ntermolecul ar 
forces and internal friction as the molecules move past each 
other. Thick fluids, like molasses, have relatively high vis- 
cosity; they do not flow easily. Thinner fluids, like water, 
have lower viscosity; they flow very easily. 

Simple kinetic theory gives a good molecular description of 
the phenomenon of viscosity in gases, but a physical description 
of the viscosity of liquids, where molecular activity is re- 
stricted, is less well developed. 

The foundations of modern viscous flow theory were laid ih 
the seventeenth century by Sir Isaac Newton, who proposed a 
method of quantifying the viscosity of a fluid. He considered 
the flow to be equivalent to a large number of thin layers of 
fluid sliding over each other. The internal friction or viscos- 
ity of the fluid was assumed to give rise to shear stresses t 
between the sliding layers. These stresses act in such a way 
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that they tend to retard the faster moving layers and accelerate 
the slower layers. Newton (1687) postulated that the viscous 
shear stresses were directly proportional to the shear strain 
rate. The relationship for internal friction in a viscous fluid 
(as proposed by Newton, 1687) can be written as 



where 

T = internal shear stress in the fluid 

n = coefficient of absolute or dynamic viscosity or coefficient 
of internal friction 
du/dz = velocity gradient across film 

It follows from equation (5.49) that the unit of viscosity 
must be the unit of shear stress divided by the unit of shear 
rate. In the newton-meter-second system, the unit of shear 
stress is the newton per square meter. Hence, the unit of vis- 
cosity will be newton per square meter multiplied by second, or 
N s/m . In the SI system the unit of pressure or stress (N/m ) 
is known as a pascal, abbreviated Pa, and it is becoming in- 
creasingly common to refer to the SI unit of viscosity as the 
pascal second (Pa s). In the cgs system, where the dyne is the 
unit of force, viscosity is expressed as dyne-second per square 
centimeter. This unit is called the poise after Poiseuille who, 
though primarily interested in the movement of olood, studied 
the flow of distilled water in very narrow tubes. ,It happens 
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that a viscosity of 1 poise is rather high; hence the frequent 
use of the centipoise (cP). 

Conversion of viscosity from. one system to another can be 
facilitated by Table 5.1. To convert from a unit in the column 
on the left side of the table to a unit at the top of the table, 
multiply by the corresponding value given in the table. For 
example, n = 0.04 N s/m^ = 0.04x1.45x10“^ Ibf s/in^ = 

5.8x10“^ Ibf s/in^. Three metric and one English system are 
presented - all based on force, length, ana time. Metric units 
are the centipoise, the kilogram-second per square meter, and 
the newton-second per square meter (or Pa s). The English unit 
is pound-second per square inch, or reyn, in honor of Osborne 
ReynolcLs. 

Viscosity is the most important property of the luoricants 
employed in hydrodynamic and elastohydrodynamic lubrication. In 
general, however, the viscosity of a lubricant does not simply 
assume a uniform level in a given bearing. This results from 
the nonuniformity of the pressure and/or temperature prevailing 
in the lubricant film. Indeed many elastohydrodynamical ly lu- 
bricated machine elements operate over ranges of pressure and/or 
temperature so extensive that the consequent variations in trie 
viscosity of the lubricant may become substantial and, in turn, 
may dominate the operating characteristics of the bearing. 

Consequently an adequate knowledge of the viscosity-pres- 
sure and viscosity-pressure-temperature relationships of lubri- 
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cants is indispensable. The next two sections therefore deal 
with such relationships. 

5.4.1 Effect of Pressure on Viscosity 

As long ago as 18y3, Barus proposed the following formula 
for the isothermal vi scosity-pressure dependence of liquids: 

(5.50) 

where 

n = viscosity at gauge pressure 
rig = viscosity at atmospheric pressure 
a = pressure-viscosity coefficient of lubricant 
The pressure-viscosity coefficient a characterizes the liquid 
considered and depends only on temperature, not on pressure. 
Although equation (5.50) is extensively used, it is not gener- 
ally applicable and is valid as a reasonable approximation only 
in a moderate-pressure range. 

Because of the shortcomings of equation (5.50), several 
isothermal viscosity-pressure formulas have been proposed that 
usually contain two or more' parameters instead of the single 
parameter suggested oy Barus (1893). One of these approaches, 
which is used in this book, was developed by Roelands (1966), 

•k 

log denotes the common or Briggsian logaritnm, log^g; 

In denotes the natural or Napierian logarithm, log^. 
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who undertook a wide-ranging study of the effect of pressure on 
the viscosity of lubricants. For isothermal conditions the 
Roelands (1966) formula can he written as 

log n + 1.200 = (log no + 1.200)(l + (b.bl) 

where 

p = gauge pressure, kgf/cm'^ 

= viscosity-pressure index, a dimensionless constant 
By taking the antilog of both sides and rearranging tenms, 
equation (5.51) becomes 

z z 

(l+p/2000) 1 1.2[(l+p/2000) ^ - ll 

n = Oq X 10 ^ 

By writing this in dimensionless terms and rearranging, 

[l-(l+p/c)^l] 

(6.52) 

where 

= 6.31x10“^ N s/m^ (9.15x10"^ lof s/in^') 
c = 1.96x10® N/m^ (28,440 Ibf/in®) 

In equations (5.51) and (5.52), care must be taken to ensure 
that the same dimensions are used in defining the constants. 

In the Roelands (1966) formulation the lubricant is defined 
by the atmospheric viscosity iiq, the viscosity-pressure in- 
dex Z^, and the asymptotic isoviscous pressure p^^ 

The equation describing the latter quantity can be written as 
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Blok (1965) arrived at the very important conclusion that all 
elastohydrodynamic lubrication results achieved hitherto for an 
exponential-pressure dependence (see equation (5.50)) can, to a 
fair approximation, be generalized simply by substituting the 
reciprocal of the asymptotic isoviscous pressure 1/P^^ 
for the viscosity-pressure coefficient a occurring in those 
results. This implies that 

a«— ^ (5.54) 

^iv,as 


5.4.2 Effects of Pressure and Temperature on Viscosity 

The viscosity is found to be extremely sensitive to both 
pressure and temperature. This extreme sensitivity forms a con- 
siderable obstacle to the analytical description of the conse- 
quent viscosity changes. Roelands (1966) noted that at constant 
pressure the viscosity increases more or less exponentially with 
the reciprocal of absolute temperature. Similarly at constant 
temperature the viscosity increases more or less exponentially 
with pressure as shown, earl ier in this chapter. In general, 
however, the relevant exponential relationships constitute only 
first approximations and may be resorted to only in moderate- 
temperature ranges. 

From Roelands (1966) the viscosity-temperature-pressure 
equation can be written as 
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log (log n + 1.200) = - Sq log [1 + 


‘■m 


135 




m 


135 


log (l + — 2 — ) + log G„ 

V 2000/ ^ 


Taking the antilog we get 


log n + 1.200 - G 


y-^ 2000/ 


-C2log(l+T^/l35)+D2 


1 + 


_TnL^^° 
135 


(5.55) 


According to equation (5.55) four parameters Gq, Sq, C 2 , and 
are sufficient to enable the viscosity n to be expressed 
in centipoise as a function of temperature in degrees 

Celsius and gauge pressure p in kgf/cn/. 


5.5 Density 

The effects of temperature on viscosity were found in the 
preceding section to be most important. For a comparable change 
in pressure and/or temperature the density change is small com- 
pared with the viscosity change. However, very high pressures 
exist in elastohydrodynamic films, and the liquid can no longer 
be considered as an incompressible medium. It is therefore 
necessary to consider the dependence of the density on pressure. 

The variation of density with pressure is roughly linear at 
low pressures, but the rate of increase falls away at high pres- 
sures. The limit of the compression of mineral oils is only 
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25 percent for a maximum density 

From Dowson and Higginson (1966) 

mineral oil can be written as 

P = — = 
PQ 


increase of about 33 percent, 
the dimensionless density for 


1 + P 
1 + 1.7 p 


(5.56) 


where 

Pq = density at atmospheric conditions 
p = gauge pressure, GPa 

Therefore the general expression for the dimensionless density 
can be written as 


p 


1 + 


APE' 

1 + BPE' 


(5.57) 


where A and B are constants dependent on the fluid. 


5.6 Energy Equation 


The equation used to determine the distribution of tempera- 
ture within the fluid is a mathematical statement of the princi- 
ple of energy conservation. As the lubricant is sheared, work 
is done on it, and there is a temperature rise that in turn 
changes the viscosity of the fluid. This variation of viscosity 
must be included in many solutions of the Reynolds equation, 
particularly when sliding occurs. Likewise from the standpoint 
of heat transfer and thermal distortion it is desirable to de- 
termine the temperature gradients of the energy equation for a 
fluid, and the full equation can be written as 
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p 


dE, 

c 

It 


. = J*' 


3x 



8y/ 



Convection Conduction 



- P 


9u 3v 9w\ 
9x 9y 3z/ 


+ Y 


Adiabatic 

compression 

Viscous 

dissipation 


(b.ba) 

where 



+ 2 


9u 9w _9w _l_ _A 

9z 9x 9x 9y 9y dz) 3 





4. .iX ( _ -^1 4. 

9y y9y 9xy 9z ^9z 

The convection term can usually be dropped from the energy 
equation in an analysis of elastohydrodynamic lubricating 
films. Likewise the rate at which heat can be conducted along 
the film (i.e., in the x and y directions) is small compared 
with the rate of conduction to the bearing surfaces. The sig- 
nificant terms in the viscous dissipation expression are seen to 
be those involving differentials of u and v in the 2 di- 
rection, and the reduced energy equation for lubricating films 
thus becomes 



(b. 60 ) 


9v 

9y. 


(b.b9) 
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If the lubricant is considered to be an incompressible 
fluid, the adiabatic compression term in equation (5.60) is 
zero. Furthermore the thermal conductivity of liquid lubricants 
can often be treated as a constant in lubrication problems, and 
if both these assumptions are acceptable, the energy equation 
adopts the following form: 


J*ic 


92t 


m 


9z^ 


+ n 



(5.61) 


This equation represents a balance between ttie rate of heat pro- 
duction within the film by viscous action and the rate of heat 
conduction from the film to the solids. 


5.7 Elasticity Equation 

The next task is to calculate the deformation in the ellip- 
soidal solids that represent an elastohydrodynamic contact as 
shown in Figure 2.18. Dowson (1965) distinguished between two 
modes of deformation that may exist in machine elements. In one 
mode the contact geometry may be affected by overall distortion 
of the elastic machine element resulting from applied loads, as 
shown in Figure 5.3(a). In the other the normal stress distri- 
bution in the vicinity of the contact zone may produce local 
elastic deformations that are significant when compared with the 
lubricant film thickness, as shown in Figure 5.3(b). This is 
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the mode of deformation with which this book is concerned. The 
important distinction is that the first form of deformation is 
relatively insensitive to the distribution and magnitude of the 
stresses in the contact zone, whereas the second mode of defor- 
mation is intimately linked to the local stress conditions. 

The correct evaluation of elastic deformation on the sur- 
face of a solid depends on an adequate representation of the 
applied normal pressures. The simplest procedure is to divide 
the actual pressure distribution into rectangular blocks of uni- 
form pressure and to permit each rectangle to be of such small 
dimensions that adequate predictions of elastic displacements 
ensue. More complex representations of two-dimensional pres- 
sure distributions within each rectangle would generally permit 
larger rectangles to be used, but the additional complexity of 
the expressions and the added computation times involved in nu- 
merical solutions make it desirable to exploit the simpler rep- 
resentation to the fullest. 

The deformation analysis is developed in general form here 
since it is used in Chapter 7 in developing the elastohydrody- 
namic lubrication tneory. The Dowson and Hamrock (1976) work is 
used extensively in this section. 

In Chapter 3 the general geometry of two ellipsoidal solids 
in elastic contact is described. In the subsequent analysis it 
is convenient to consider the deformation of an equivalent elas- 
tic half-space subjected to a pressure distribution over the 
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ellipse of semimajor and semiminor axes a and b as previous- 
ly defined in equations (3.13) and (3.14). The resulting 
elastic deformation can be considered to be equivalent to the 
total deformation of two elastic ellipsoids having elastic 
constants and respectively, if the half-space is 

allocated the equivalent elastic parameter E' defined by equa- 
tion (3.16). 

Figure 5.4 shows a rectangular area of uniform pressure 
with the coordinate system to be used. From Timoshenko and 
Goodier (1951) the elastic deformation at a point (x,y) of a 
semi-infinite solid subjected to a pressure p at the point 
(^l,yi) can be written as 

2p dx 

d6 = 

irE'r 

The elastic deformation at a point ("x,y) due to the uniform 
pressure over the rectangular area 2a x 2b is thus 


a b 



Integrating the preceding equation gives 

5 = f PD (5.62) 

where 
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D = (x + b) In 


_(y + a) + 1 

-(y + a)^ + Cx + bid^^^ 

(y - a) + 

[(y - a)^ + (x + b)^_ 

1/2 


+ Cy + a) In 


+ (x - b) In 


+ (y - a) In 


(5E + ¥) + 

-Cy + + Cx + b)^ 

1/2 

_ (x - b) + 

[cy + a)^ + (x - b)^ 

1 
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- -2 _ -2 
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(y + "a) + 

Cy + a)^ + Cx - b)2 

[ 1/2 

(x - b) + 

- 2 - - ^ 
.Cy - a) + Cx - b)J 

1/2^ 

_ (x + b^) + 

Cy - a)^ + Cx + F)^ 

1/2 


(5.63) 

As a check on the validity of equation (5.62) the following 
two cases were evaluated: 

Case 1 : For b = a and x = y = 0, equation (5.62) reduces to 

/ 

6 = ^ pI ln(l + */^) ( 6 . 64 ^ 

Equation (5.64) represents the elastic deformation at the center 
of the square of uniform pressure. This equation is in exact 
agreement with that presented by Timoshenko and Goodier (iybi). 
£ase_2: For b = a and x = y = a, equation (5.62) reduces to 
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5 


!pI 


ln(l + 




(6.bb) 


Equation (5.65) represents the elastic deformation at the corner 
of a square of uniform pressure. This equation is also in 


agreement with the result presented by Timoshenko and Goodier 
(1951). From equations (5.64) and (5.65) we find the corner 
deformation to be one-half the deformation at the center of a 
square block of uniform pressure. 

Now the term & in equation (5.62) represents the elastic 
deformation at a point (x,y) due to a rectangular area 2a x 2^ 
of uniform pressure p. If the contact ellipse is divided into 
a number of equal rectangular areas, the total deformation at a 
point (x,y) due to the contributions of the various rectangular 
areas of uniform pressure in the contact ellipse can be evalua- 
ted numerically. Figure 5.5 shows how the area inside and out- 
side the contact ellipse can be divided into a number of equal 
rectangular areas. For illustration the contact is shown divi- 
ded into a grid of 6 x 6 rectangular areas. The arrangement 
illustrated in Figure 5.5 can be used to evaluate the total 
elastic deformation, caused by the rectangular areas of uniform 
pressure within the contact ellipse, at any point inside or out- 
side the contact ellipse as 





6 ~ 



1 = 1 , 2 ,.. 


^1, j^m,s 


(5.66) 
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where 


m = |k - i| + 1 


(5.b7) 


s = - j| + 1 (5.68) 

Note that ^ would be 0 in equation (5.63) evaluated at 

X = 0, y = 0, while O 2 3 would be evaluated at x = 2b, 7 = 47. 

The elastic deformation at the center of the rectangular 
area 6g g shown in Figure 5.5, caused by the pressure of the 
various rectangular areas in the contact ellipse, can be written 
as 


'^9.5 = 7T ri.1^9.5 + ^2,1°8,5 




^ 1 , 2 ® 9,4 ^ 2 , 2 ° 8,4 + '•* + ^ 6 , 2 ® 4,4 


^ 1 , 6 ^ 9, 2 ^2 


.Dc 


+ 


+ P 




(5.69) 


Equation (5.69) points out more explicitly the meaning of 
equation (5.66). A numerical analysis of the elastic deforma- 
tion of a contacting ellipsoidal solid and plane was performed 
by Dowson and Hamrock (1976). The analysis assumed that the 
pressure in the contact zone was Hertzian. It also assumed that 
the contact zone could be divided into rectangular areas witn 
uniform pressure within each rectangular area. The resulting 
equations were programmed on a digital computer. Four limiting 
conditions were evaluated. They consisted of two extremes of 
applied normal load: a light load of 8.964 N (2 Ibf), and a 
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heavy load of 896.4 N (200 Ibf). The two other extremes are of 
the curvature of the contacting solids: two equal spheres in 

contact, and a ball and outer race of a ball bearing. It was 
speculated that conclusions drawn from the results obtained for 
these limiting conditions could also be made for any interme- 
diate condition. 

Figures 5.6, 5.7, and 5.8 illustrate the results of an in- 
vestigation of the effects of mesh size on the accuracy of the 
calculations of elastic deformation presented by Dowson and 
Hamrock (1976). In these figures the solid curve represents the 
case of equal spheres in contact in which = 0.005558 m 

(0.2188 in.), and the dashed curve represents the contact be- 
tween a ball and outer race in a bearing in which R = 0.0128 m 
(0.5055 in. ) and by R^ = 0.150 m (5.91 in. ). 

Figure 5.6 shows how the percentage difference in elastic 
deformation varies along the semimajor axis when m = 3, 4, 
and 5 - the predictions being related to more exact solutions 
for three times as many divisions of the axes in each case 
(i.e., m = 9, 12, and 15, respectively). In this figure, we see 
a large drop in the error term 100 [(6- - 63i^)/63-J when 
going from m = 4 to m = 5. It was therefore argued that a 
value of trT of 5 would yield acceptable accuracy with minimum 
computing effort under most conditions. 

The ratio of the elastic deformation to the distance separ- 
ating the two undeformed solids in contact is shown as a func- 
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tion of the distance along the semimajor and semiminor axes in 
Figures 5.7 and 5.8, respectively. These figures show that the 
elastic deformation ceases to be a significant component of the 
overall separation of the solids at modest distances away from 
the edge of the Hertzian contact zone. To be more specific, 
from the curves shown in Figures 5.7 and 5.8 we see that for 
equal spheres in contact, represented by solid lines in the 
figures, 6/S < 0.05 corresponds to x > 2.6 b and y > 2.6 a. 
This means that the elastic deformation is less than 5 percent 
of the film thickness due to the geometry of the undeformed 
solids at distances of only 1.6 times the semiaxes dimensions 
beyond the edge of the Hertzian contact zone. For the ball-and- 
outer-race example a value of 6/S < 0.05 is achieved when 
y > 1.9 a and x > 4.0 b. This means that the elastic deforma- 
tion is less than 5 percent of the film shape due to the unde- 
formed geometry at distances beyond the contours of the Hertzian 
ellipse of only 0.9 and 3.0 times the semimajor and semiminor 
axes, respectively. 

The great simplicity introduced into the calculation of 
elastic deformations by the assumption that the pressure is uni- 
form over rectangular blocks has already been mentioned. This 
assumption enabled the wide range of approximate numerical solu- 
tions required to generate the empirical film thickness equa- 
tions to be obtained in an acceptable time and with a justifi- 
able computing effort. Evans and Snidle (1978) have discussed a 
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refinement to the numerical procedure based largely on a more 
detailed approach to the calculation of elastic deformations. 

In this approach tne elastic defonmation at any point resulting 
from a block of pressure in region A is determined by dividing A 
into two regions, B and C, and adopting different procedures for 
the evaluation of the basic integrals in the two regions. In 
the central region, C, where a singularity occurs in a straight- 
forward approach to the problem, the pressure function is ap- 
proximated by a biquadratic polynomial to facilitate numerical 
analysis. In the remaining region, B, the integral is evaluated 
in a straightforward manner by the application of Simpson's rule. 

Evans and Snidle (1978) provided two full solutions to the 
elastohydrodynamic problem for nominal point or circular con- 
tacts (k = 1) that revealed most clearly all the characteristics 
of film shape and pressure distribution associated with the 
problem. They also investigated the influence of lubricant 
starvation, as determined by the location of the inlet boundary 
and the location of the side boundary to the computing zone, on 
film thickness. They concluded that both factors could play an 
important part in determining the accuracy of numerical solu- 
tions, and both factors are considered in detail in the numeri- 
cal solutions presented in this text. 

/ 

A comparison between the predictions of the film thickness 
equations developed in this text, which are based on solutions 
of the elasticity equation as outlined in this chapter, and the 
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refined approach presented for only two conditions by Evans and 
Snidle (1978) provides justifiable confidence in the adequacy of 
the numerical procedures adopted. In one particular case, cor- 
reponding to a sphere of radius 0.025 m rolling near a plane 
with a velocity of 0.5 m/s while supporting a load of 120 N, the 
central and minimum film thicknesses calculated by Evans and 
Snidle (1978) were predicted to within 0.3 and 1.4 percent, re- 
spectively, by equations (8.41) and (8.23) in this text. We are 
grateful to Dr. Snidle (private communication, April 1980) for 
the information that enabled this comparison to be made. 

5.8 Closure 

This chapter has dealt with the basic lubrication and elas- 
ticity equations used in the study of hydrodynamic lubrication. 
The Navier-Stokes equations have been derived in general form 
and then reduced within the framework of assumptions that are 
deemed to be valid for lubricating films. The Reynolds equation 
was then derived in two different ways: one from the coupling 

of the reduced Navier-Stokes equations with the continuity equa- 
tion, the other from first principles and the direct combination 
of well-established viscous flow expressions with a statement of 
flow continuity. The Reynolds equation contains viscosity and 
density as parameters. These properties depend on the tempera- 
ture and pressure, and it is therefore necessary to establish 
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expressions relating these quantities. The energy equation, 
which governs the temperature distribution within the lubricant, 
is also presented in this chapter. It is noted that solutions 
to both the energy and Reynolds equations may be required in 
some situations. 

The Reynolds equation also contains the film thickness as a 
parameter. The film thickness is a function of the undeformed 
geometry and of the elastic behavior of the bearing surfaces. 

The elasticity equation has therefore been developed in this 
chapter. The coupling of the Reynolds and elasticity equations 
will be dealt with in Chapter 7, where the full theory of elas- 
tohydrodynamic lubrication is developed. 
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SYMBOLS 
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A*, B*, C*, 
D*, L*, M* , 
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c 


D 

D 

D 

De 

d 

"d 


d 

a 




d 


0 


constant used in equation (3.113) 
relaxation coefficients 

drag area of ball, 

semimajor axis of contact ellipse, m 
a/2m 

total conformity of bearing 
semiminor axis of contact ellipse, m 
b/2m 

dynamic load capacity, N 
drag coefficient 
constants 

19,609 N/cm^ (28,440 Ibf/in^) 

number of equal divisions of semimajor axis 

distance between race curvature centers, m 

material factor 

defined by equation (5.63) 

Deborah number 
ball diameter, m 

number of divisions in semiminor axis 
overall diameter of bearing (Figure 2.13), m 
bore diameter, m 
pitch dameter, m 

pitch d’-ameter after dynamic effects have acted on ball, m 
inner-race diameter, m 
outer-race diameter, m 
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E 


E' 


E 

E 

S 

~i 

e 

F 


F* 

F 

? 


max 


Sr 

f 


G* 

G 

g 


modulus of elasticity, N/m 


2 


effective elastic modulus, 




7 7 

internal energy, m /s 
processing factor 

^^^min “ ^min^^^min^ ^ 

2 

elliptic integral of second kind with modulus (1 - 1/k ) 

approximate elliptic integral of second kind 

dispersion exponent 

normal applied load, N 

normal applied load per unit length, N/m 

lubrication factor 

integrated normal applied load, N 


1/2 


centrifugal force, N 

maximum normal applied load (at =0), N 

applied radial load, N 

applied thrust load, N 

normal applied load at angle i}), N 

P 1/2 

elliptic integral of first kind with modulus (1 - l/k*^) 

approximate elliptic integral of first kind 

race conformity ratio 

rms surface finish of ball, m 

rms surface finish of race, m 

dimensionless materials parameter, aE 

p 

fluid shear modulus, N/m 
hardness factor 

2 

gravitational constant, m/s 
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9v 

H 

A 

H 

H_ 


c,s 


mm 


min,r 


H 


min,s 


»min 


«min 


h 

h. 


h. 


8/3 2 

dimensionless elasticity parameter, W ' /U 

3 2 

dimensionless viscosity parameter, GW /U 
dimensionless film thickness, h/R^ 
dimensionless film thickness, H(W/U)^ = F^h/u^nQR^ 
dimensionless central film thickness, h^/R^ 
dimensionless central film thickness for starved 
lubrication condition 
frictional heat, N m/s 

dimensionless minimum film thickness obtained from EHL 
elliptical-contact theory 

dimensionless minimum film thickness for a rectangular 
contact 

dimensionless minimum film thickness for starved 
lubrication condition 

dimensionless central film thickness obtained from . 
least-squares fit of data 

dimensionless minimum film thickness obtained from 
least-squares fit of data 

dimensionless central-film-thickness - speed parameter, 

h^u-o*5 

dimensionless minimum-film-thickness - speed parameter, 

new estimate of constant in film thickness equation 
film thickness, m 
central film thickness, m 
inlet film thickness, m 
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h . 


rmn 


J* 

J 

K 

k 

k 

k 



L 


50 


l 

T 

M 


film thickness at point of maximum pressure, where 
dp/dx = 0, m 

minimum film thickness, m 
constant, m 

diametral interference, m 

2 

ball mass moment of inertia, m N s 
integral defined by equation (3.76) 
integral defined by equation (3.75) 
function of k defined by equation (3.8) 
mechanical equivalent of heat 

2 

polar moment of inertia, m N s 
load-deflection constant 
ellipticity parameter, a/b 
approximate ellipticity parameter 
thermal conductivity, N/s 'C 
lubricant thermal conductivity, N/s °C 
fatigue life 
adjusted fatigue life 

reduced hydrodynamic lift, from equation (6.21) 
lengths defined in Figure 3.11, m 

fatigue life where 90 percent of bearing population will 
endure 

fatigue life where 50 percent of bearing population will 
endure 

bearing length, m 

constant used to determine width of side-leakage region 
moment, Nm 
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g 

p 
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m* 


m 
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N 

n 

n* 

P 

Po 


Hz 


max 


iv,as 


m 


gyroscopic moment, Nm 

-0 75 

dimensionless load-speed parameter, WU 
torque required to produce spin, N m 
mass of Dali, N s /m 

dimensionless inlet distance at boundary between fully 
flooded and starved conditions 
dimensionless inlet distance (Figures 7.1 and 9.1) 
number of divisions of semimajor or semiminor axis 
dimensionless inlet distance boundary as obtained from 
Wedeven, et al. (1971) 
rotational speed, rpm 
number of balls 
refractive index 

constant used to determine lengtn of outlet region 

dimensionless pressure 

dimensionless pressure difference 

diametral clearance, m 

free endplay, m 

2 

dimensionless Hertzian pressure, N/m 
2 

pressure, N/m 

2 

maximum pressure within contact, SF/Znab, N/m 

2 

isoviscous asymptotic pressure, N/m 
solution to homogeneous Reynolds equation 
thermal loading parameter 

2 

dimensionless mass flow rate per unit width, qtiQ/PQE'R 
reduced pressure parameter 

volume flow rate per unit widtn in x direction, f/ls 
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2 

volume flow rate per unit width in y direction, m /s 
curvature sum, m 

arithmetical mean deviation defined in equation (4.1), 
operational hardness of bearing material 
effective radius in x direction, m 
effective radius in y direction, m 
race curvature radius, m 

radii of curvature, m 

cylindrical polar coordinates 

spherical polar coordinates 

defined in Figure 5.4 

geometric separation, m 

geometric separation for line contact, m 

empirical constant 

shoulder height, m 

^O^^max 

tangential (traction) force, N 
temperature, “C 
ball surface temperature, “C 
average lubricant temperature, *C 
ball surface temperature rise, °C 

viscous drag force, N 
time, 5 

auxiliary parameter 

velocity of ball-race contact, m/s 
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u 

u 

TT 

AU 

V 

W 


w 

X 

Y 



a 


B 

6 ' 


r 


Y 

Y 


velocity of ball center, m/s 

dimensionless speed parameter, tiqU/E'R^ 

surface velocity in direction of motion, (u^ + Ujj)/2, m/s 

number of stress cycles per revolution 

sliding velocity, u^^ - U|^, m/s 

surface velocity in transverse direction, m/s 

2 

dimensionless load parameter, F/E'R 

surface velocity in direction of film, m/s 

dimensionless coordinate, x/R^^ 

dimensionless coordinate, y/R^ 

dimensionless grouping from equation (6.14) 

external forces, N 

constant defined by equation (3.48) 

viscosity pressure index, a dimensionless constant 

coordinate system 

• 2 

pressure-viscosity coefficient of lubrication, m /N 

radius ratio, Ry/R^ 

contact angle, rad 

free or initial contact angle, rad 

iterated value of contact angle, rad 

curvature difference 

viscous dissipation, N/m s 

total strain rate, s“^ 

elastic strain rate, s“^ 

viscous strain rate, s"^ 
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6 

6 * 

6n 


'^2 

c 



0 

A 

X 


flow angle, deg 

total elastic deformation, m 

0--1 

lubricant viscosity temperature coefficient, C 
elastic deformation due to pressure difference, m 
radial displacement, m 
axial displacement, m 
displacement at some location x, m 
approximate elastic deformation, m 
elastic deformation of rectangular area, m 
coefficient of determination 
strain in axial direction 
strain in transverse direction 
angle between ball rotational axis and bearing 
centerline (Figure 3.10) 
probability of survival 

2 

absolute viscosity at gauge pressure, N s/m 
dimensionless viscosity, ti/hq 

2 

viscosity at atmospheric pressure, N s/m 

6.31x10“^ N s/m^(0.0631 cP) 

angle used to define shoulder height 

film parameter (ratio of film thickness to composite 

surface roughness) 

equals 1 for outer-race control ana 0 for inner-race 
control 

second coefficient of viscosity 

Archard-Cowking side-leakage factor, (1 + 2/3 a^)~^ 
relaxation factor 
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p 

V 

5 

p 

"p 


^0 


"1 



T 


e 


T 


L 




$1 

ij' 

'I'o 


fir 


coefficient of sliding friction 
p/n 

Poisson’s ratio 

divergence of velocity vector, (3u/3x) + (av/ay) + (aw/az), s ^ 

2 , 4 

lubricant density, N s /m 
dimensionless density, p/pq 

2 4 

density at atmospheric pressure, N s /m 
2 

normal stress, N/m 

2 

stress in axial direction, N/m 

2 

shear stress, N/m 

2 

maximum subsurface shear stress, N/m 
2 

shear stress, N/m 

2 

equivalent stress, N/m 

2 

limiting shear stress, N/m 

ratio of depth of maximum shear stress to semiminor axis of 

contact ellipse 
Ph3/2 

auxiliary angle 
tnermal reduction factor 
angular location 
limiting value of '!> 

absolute angular velocity of inner race, rad/s 
absolute angular velocity of outer race, rad/s 
angular velocity, rad/s 

angular velocity of ball-race contact, rad/s 
angular velocity of ball about its own center, rad/s 
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Subscripts: 

a 

b 


c 

be 

IE 

IR 

i 

K 

min 

n 

0 

PVE 

PVR 

r 

s 

x,y,z 

Superscript: 

(~) 


angular velocity of ball around shaft center, rad/s 
ball spin rotational velocity, rad/s 

solid a 
solid b 
central 
ball center 

isoviscous-elastic regime 

isovi scous-rigid regime 

inner race 

Kapitza 

minimum 

iteration 

outer race 

piezovi scous-elastic regime 
piezoviscous-rigid regime 
for rectangular area 
for starved conditions 
coordinate system 

approximate 
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Figure 5. 5. - Division of area in and 
around contact into equal rectangular 
areas. 
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Figure 5. 7. - Effect of location along semimajor axis on ratio of 
elastic deformation to distance se^rating two solids in contact 
due to geometry of the solids for m • 5. 





Figure 5,8. - Effect of location along semiminor axis on ratio 
of elastic deformation to distance separating two solids in 
contact due to geometry of the solids for m • 5. 
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